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摘 要 : 中 国 西北 干旱 区 是 典型 的 地 下 水 依赖 型 生态 系统 ,地 下 水 资源 是 西北 干旱 区 发 展 的 战略 性 支撑 ,实时 掌 
尾 地 下 水 的 动态 变化 对 揭示 地 下 水 的 变化 特征 以 及 区 域 经 济 发 展 具有 重要 意义 。 以 往 基 于 人 工 的 传统 地 下 水 鉴 
测 研究 耗费 巨大 ,难以 实现 大 尺度 长 时 间 序 列 的 分 析 , 重 力 卫 星 的 运用 为 解决 这 一 问题 提供 了 途径 ,过 往 研 究 尚 未 
对 西北 干旱 区 地 下 水 资源 量 可 持续 性 进行 量化 评价 ,为 探究 近年 来 西北 干旱 区 地 下 水 资源 量 可 持续 性 情况 ,基于 
GRACE 重力 卫星 和 GLDAS 数据 反 演 中 国 西北 干旱 地 区 2002—04—2020—09 共计 200 个 月 的 地 下 水 储量 变化 情 
bL ,判断 其 趋势 变化 ,再 计算 基于 GRACE 的 地 下 水 干旱 指数 获取 西北 干旱 区 地 下 水 干旱 情况 ,最 后 进一步 将 GGDI 
与 可 持续 性 指数 相 结合 ,量化 地 下 水 可 靠 性 . 回 弹性 脆弱 性 ,以 此 评价 近年 来 中 国 西北 干旱 区 地 下 水 资源 量 可 持 
续 性 。 结 果 显 示 :2002 一 2020 年 ,西北 干旱 区 地 下 水 储量 大 幅度 减少 ,地 下 水 干旱 程度 加 深 , 地 下 水 平均 可 靠 性 为 
0.495 , 回 弹性 0.470 ,脆弱 性 为 0.404, 区 域 的 地 下 水 资源 量 可 持续 性 指数 SI 为 0.28 ,属于 较 低 的 地 下 水 资源 量 可 持 
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续 性 水 平 。 研 究 揭 示 了 西北 干旱 区 的 地 下 水 资源 量 可 持续 性 时 空 分 布 规律 及 演变 趋势 ,为 区 域 地 下 水 资源 保护 和 


管理 提供 理论 和 数据 支撑 。 


关键 词 : 西北 干旱 区 ; 重力 卫星 ; 地 下 水 干旱 指数 ;地 下 水 资源 量 ; 


地 下 水 作为 水 资源 的 重要 组 成 部 分 ,具有 分 布 
范围 广 水 质 优良 、 开 发 成 本 低 等 特点 "3。 中 华人 
民 共 和 国 水 利 部 发 布 的 2020 年 度 《 中 国 水 资源 公 
报 》 指 出 ,全 国 地 下 水 资源 总 量 为 8553.5x10? m, d 
水 资源 总 量 的 27% ,2020 年 地 下 水 资源 供给 量 为 
892.5x10* m, 占 供水 总 量 的 15.4%, 地 下 水 资源 在 供 
给 中 占有 极其 重要 的 地 位 “5 。 随 着 用 水 需求 的 不 
断 提升 ,21 世纪 以 来 ,我 国 地 下 水 开采 量 持 续 超 过 
1000x10* m, 地 下 水 的 超 采 势必 造成 地 下 水 储量 减 
> 地 下 水 水 位 的 下 降 中 ,可 能 导致 地 下 水 依赖 型 生 
态 系统 退化 并 危及 区 域 生 态 安全 。 干 旱 半 干旱 区 
就 是 典型 的 地 下 水 依赖 型 生态 系统 ,只 有 可 持续 
利用 地 下 水 ,才能 缓解 干旱 半 干 旱 区 地 下 水 环境 
的 恶化 "1。 

地 下 水 可 持续 性 的 概念 最 早 由 Alley 等 中 在 
1999 年 提出 ,后 经 过 不 断 的 完善 可 以 定义 为 :地 下 
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可 持续 性 


水 系统 能 够 为 生态 系统 和 人 类 社会 持续 提供 足够 
数量 和 良好 质量 的 水 资源 的 能 力 。 国 内 外 也 展开 
了 大 量 的 地 下 水 可 持续 性 相关 的 研究 ,主要 包括 基 
于 数量 的 、 基 于 质量 的 以 及 综合 评价 方法 3 种 。 在 
地 下 水 数量 评价 方面 ,Thomas 等 "通过 GRACE 
(Gravity Recovery and Climate Experiment) 和 GLDAS 
(Global Land Date Assimilation System ) 水 文 模型 , 结 
合 基于 GRACE 的 地 下 水 干旱 指数 (CRACE Ground- 
water Drought Index, GGDI) , 评 价 了 2002 一 04 一 
2016 一 07 全 球 37 个 大 型 含水 层 的 地 下 水 资源 量 可 
持续 性 。 地 下 水 质量 评价 方面 主要 包括 通过 地 下 
水 指数 .地 下 水 脆弱 性 以 及 地 下 水 质 可 持续 性 指标 
进行 评价 。 综 合 评 价 是 指 综合 考虑 地 下 水 的 数量 
和 质量 来 评价 地 下 水 可 持续 性 ,如 Singh 等 "通过 
构建 基于 15 个 指标 的 地 下 水 可 持续 性 指数 GSI 评 
估 印 度 拉 贾 斯 坦 邦 西部 干旱 地 区 地 下 水 可 持续 性 ; 
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Jiang 等 由 在 驱动 力 -- 压 力 -- 状 态 -- 影 响 -- 响 应 框 
架 的 基础 上 ,建立 了 地 下 水 可 持续 性 指标 体系 ,并 
运用 层次 分 析 法 和 信息 粹 理论 重新 构建 了 地 下 水 
可 持续 性 评价 方法 。 

开展 地 下 水 的 实时 动态 监测 是 地 下 水 资源 量 
可 持续 性 利用 的 重要 保障 ,传统 的 地 下 水 储量 监测 
技术 主要 是 通过 监测 井 来 实测 地 下 水 位 数据 以 达 
到 估算 地 下 水 储量 的 目的 ,然而 该 方法 存在 较 大 局 
限 ,包括 水 井 的 数量 少 、 分 布 不 均 以 及 很 难 对 大 尺 
度 的 区 域 进行 监 测 等 汪 ”。 随 着 重力 卫星 探测 技术 
的 不 断 发 展 , 越 来 越 多 的 国内 外 研究 者 通过 利用 
GRACE 卫星 数据 ,更 好 的 实现 对 大 尺度 区 域 地 下 水 
的 监测 。Neves 等 中 评估 比较 了 几 种 不 同 来 源 的 
GRACE 数据 在 伊比 利 亚 地 区 进行 地 下 水 储量 监测 
中 的 性 能 ,结果 显示 CSR Mascon 数据 在 该 地 区 表现 
最 好 ;Karunakalage 等 5 利用 降低 尺度 的 高 分 辩 率 
GRACE 数据 监测 印度 古吉拉特邦 及 周边 地 区 的 地 
下 水 储量 减少 情况 并 验证 了 GRACE SH 数据 在 水 管 
理 人 研究 中 的 适用 性 ; 束 秋 妍 等 “利用 GRACE 重力 
卫星 数据 反 演 华北 平原 地 下 水 储量 变化 探索 华北 
平原 地 下 水 储量 变化 的 时 空 特征 和 驱动 机 制 ; 孙 依 
SF EP GRACE 卫星 结合 GLDAS 水 文 模型 对 和 田 
地 区 克 里 雅 河流 域 11 a 间 的 陆地 水 储量 动态 变化 进 
行 反 演 ,模拟 地 下 水 等 效 水 柱 高 变化 趋势 以 此 构建 
地 下 水 水 位 估算 模型 ,结合 实测 数据 验证 精度 , 验 
证 结果 表明 模型 适用 性 较 好 ; 陶 征 广 等 所 基于 
CRACE 和 GLDAS 反 演 安 徽 省 地 下 水 储量 的 时 间 变 
化 序列 , 反 演 结果 与 国家 统计 局 数据 相关 系数 达到 
89.62%。 综 上 所 述 ,GRACE 卫星 数据 在 反 演 大 尺度 
范围 地 下 水 储量 变化 有 较 高 的 精度 ,可 以 较为 系 
统 .准确 的 获取 地 下 水 储量 变化 ;并 且 GRACE 卫星 
数据 在 中 国 地 区 地 下 水 储量 监测 方面 效果 良好 。 

中 国 西北 干旱 区 降水 稀少 .蒸发 强烈 ,水 资源 
紧缺 ,自然 湖泊 湿地 天然 植 被 对 地 下 水 具有 强烈 
依赖 性 。 近 些 年 来 ,不 够 合理 的 地 下 水 资源 开发 导 
致 地 下 水 位 下 降 、 议 漠 化 入 侵 , 地 下 水 长 期 处 于 超 
采 状 态 ,地 下 水 资源 受到 严重 威胁 "1。 目 前 ,针对 
西北 干旱 地 区 地 下 水 的 研究 主要 集中 在 对 地 下 水 
储量 变化 和 地 下 水 埋 深 分 析 等 方面 ,在 地 下 水 资源 
量 可 持续 性 评价 方面 的 研究 鲜 有 报导 。 地 下 水 作 
为 西北 干旱 区 经 济 发 展 极 为 重要 的 资源 ,开展 区 域 
整体 地 下 水 可 持续 性 水 平 的 评估 ,对 区 域 发 展 、 地 


下 水 可 持续 利用 具有 重要 意义 。 本 文 旨 在 通过 利 
用 GRACE 数据 和 GLDAS 水 文 模型 反 演 西北 干旱 区 
地 下 水 储量 变化 ,并 结合 基于 GRACE 的 地 下 水 干 
旱 指 数 GGDI 进 一 步 量 化 地 下 水 资源 量 的 可 靠 性 、 
回 弹 效 .脆弱 性 ,以 此 评价 西北 干旱 区 地 下 水 资源 
量 可 持续 性 指数 ,揭示 西北 干旱 区 地 下 水 的 可 持续 
性 水 平 , 为 西北 干旱 区 地 下 水 资源 量 可 持续 性 管理 
提供 重要 参考 。 


1 研究 区 概况 与 方法 


1.1 研究 区 概况 

中 国 西北 干旱 区 位 于 亚洲 中 部 史 ,地 理 位 置 位 
T 73°~107°E, 35°~5S0°N ,东部 以 贺兰山 为 界 , 南 至 
昆仑 山 、. 阿 尔 金山 与 祁连山 , 北 侧 与 西 侧 以 国境 线 
为 界 ( 图 1)。 行政 区 上 包括 新 疆 全 境 , 甘 肃 河西 五 
市 及 内 蒙古 阿拉 善 盟 的 广大 地 区 。 区 域 总 面积 》 
2.09x10* km , 约 占 中 国 总 面积 的 5。 人 研究 区 深 居 
内 陆 ,区 域内 除了 阿勒泰 市 东北 部 , 塔 城 地 区 西北 
部 , 博 乐 市 ,奎屯 市 ` 沙 湾 县 以 及 玛 纳 斯 县 南部 , 乌 
鲁 木 齐 市 ,甘肃 河西 五 市 在 内 的 部 分 地 区 年 降水 量 
大 于 200 mm 外 ,其 他 区 域 丝 为 年 降水 量 小 于 200 
mm 的 极端 干旱 区 , 属 典型 的 大 陆 性 气候 ,地 形 以 高 
原 、 山 地 分 地 为 主 ,地 势 起 伏 大 。 其 中 ,包括 内 蒙 
二 高原. 阿尔泰 出 山脉 ,天山 山脉 .昆仑 山脉 .阿尔 
金山 脉 ,祁连山 山脉 .准噶尔 盆地 . 叶 鲁 番 盆 地 和 塔 
里 木 盆地 。 
1.2 数据 与 处 理 
1.2.1 GRACE 数据 2002 年 3 月 发 射 的 重力 恢复 
与 气候 试验 卫星 (Gravity Recovery and Climate Ex- 
periment, GRACE ) 对 全 球 重 力 场 进行 长 达 16 a 的 观 
测 ,2018 年 5 月 在 GRACE 卫星 停止 工作 1 a 之 后 
GRACE-FO (GRACE Follow-on) 发 射 成 功 开始 延续 
GRACE 卫星 的 任务 1。 

本 文选 用 德国 喷气 推进 实验 室 (JPL) 发 布 的 空 
[B] 4) 3E Fy 0.5? x 0.5? ,时 间 分 辨 率 为 1 个 月 的 
RL06-level2 Mascon 数据 模型 2 和 国家 青藏 高 原 科 
学 数据 中 心 (TPDC) 在 2020 年 10 月 发 布 的 中 国 区 域 
基于 降水 重 构 陆 地 水 储量 变化 数据 集 的 JPL RL06 
Mascon 部 分 数据 ”-, 所 有 数据 都 扣除 了 2004 一 2009 
年 间 的 平均 值 , 代 表 每 个 月 的 数据 是 相对 于 2004 一 
01 一 2009 一 12 平 均值 的 质量 变化 ,换算 为 等 效 水 柱 
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注 : 底 图 采用 自然 资源 部 标准 地 图 制作 , 审 图 号 为 GS(2016)1569 ,对 底 图 边界 无 修改 。 下 同 。 
图 1 中 国 西北 干旱 区 


Fig. 1 Arid area in Northwest China 


高 , 即 地 下 水 储量 变化 的 质量 换算 到 假想 平面 内 水 
柱 的 高 度 。 

人 研究 选择 的 时 间 跨 度 为 2002 一 04 一 2020 一 09， 

共计 200 个 月 ,部 分 月 份 缺失 ,其 中 2017 一 07 一 
2018—05 Grace 数据 间断 期 的 11 个 月 使 用 中 国 区 
域 基于 降水 重 构 陆 地 水 储量 变化 数据 集 , 其余 189 
个 月 使 用 JPL 提 供 的 GRACE/GRACE-FO 数据 。 
1.2.2 GLDAS 数据 全 球 陆 面 同 化 系统 (Global 
Land Date Assimilation System, GLDAS) 是 由 美国 宇 
航 局 哥 达 德 飞行 中 心 和 美国 国家 环境 预报 中 心 共 
同 提供 。GLDAS 数 据 包含 3 个 陆 面 过 程 模型 (CLM、 
NOAH Mosaic) ,1 个 水 文 模型 (VIC)。 数 据 包 括 了 
全 球 1979 年 至 今 的 陆 面 地 表 径 流 气温、 冠 层 水 ZK 
散发 .土壤 湿度 . 雪 水 等 数据 ,由 4 套 模型 进行 模 
40. , zs [RIA HESS Jy 19x19 5 0.259x0.25? ,时 间 分 辩 率 
为 3h 和 1 个 月 等 。 

选用 GLDAS 全球 陆 面 同化 系统 中 的 NOAH 陆 
面 过 程 模 型 来 获取 积 雪 、 冠 层 水 及 土壤 水 组 成 的 地 
表 水 储量 变化 ,时 间 分 辨 率 为 30 d, 空间 分 辩 率 为 
0.25°x0.25° 
12.3 降雨 数据 ”降雨 数据 来 源 于 国家 气象 科学 数 
据 共 享 服务 平台 中 国 地 面 气候 资料 日 值 数据 


集 (V3.0) ,包含 了 2000 一 2019 年 中 国 各 省 .市 .地 区 
的 年 月 .逐日 的 平均 降水 量 数据 (20:00 一 20:00 累 
计 降 雨量 ,单位 为 0.1 mm). 
1.2.4 土地 利用 数据 Globeland 30 是 国家 基础 地 
理 信息 中 心 研 制 的 高 分 辨 率 地表 和 覆盖 产品 ,空间 分 
JE 30 m, 分 类 精度 为 83% ,数据 在 质量 和 空间 
上 具有 较 强 的 一 致 性 ,包括 耕地 、 林 地 .草地 湿地 、 
TKS FER 灌木 地 、 人 造 地 表 ( 城 镇 等 各 类 居民 用 
地 ). 裸 地 ,冰川 和 永久 积 雪 共 10 类 。 

选用 2000 年 .2010 年 .2020 年 3 期 Clobeland 30 
数据 通过 拼接 EIER FDS HR AB IS Pr PH A F 
旱 区 土地 利用 变化 对 地 下 水 资源 量变 化 的 影响 。 
1.3 方法 
13.1 地 下 水 储量 估算 ”通过 GRACE 数据 产品 获 
得 陆地 储 水 月 变化 ATWS ,其 中 包括 了 地 上 部 分 ( 即 
植物 冠 层 表面 水 .土壤 水 和 雪 水 ) 和 地 下 部 分 ( 即 地 
下 水 ); 地 上 部 分 需要 通过 GLDAS 数 据 获取 。 

ATWS=AW. +AW.,+AW., tAGWS (1) 

SUP: AW, 为 冠 层 表 面 水 变化 ; AW a 为 土壤 水 变 
化 ; AW,,,,, 为 雪 水 变化 ; AGWS 为 地 下 水 变化 。 结 合 
GRACE 数据 和 GLDAS 数据 来 获取 地 下 水 储量 
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AGWS opace JARU P: 
AGWS crace = ATWS - (AW, AW a + AW.,,,) (2) 
1.3.2 基于 GRACE 的 地 下 水 干旱 指数 GGDI GGDI 
(Grace Groundwater Drought Index) 是 一 种 基于 归 一 
化 的 地 下 水 存储 指标 ,用 于 评估 含水 层 系统 的 地 下 
水 干旱 。 其 新 疾 性 是 其 通过 结合 地 下 水 存储 的 赤 
字 和 三 余 来 实现 对 地 下 水 干旱 的 评估 , 当 GGDI>0 
表示 盘 余 状态 ,反之 则 为 赤字 ”1。 
首先 计算 每 月 气候 (C) ,这 里 的 气候 不 是 气候 
学 定义 中 的 气候 ,气候 (C) 是 用 来 消除 月 度 变 化 因 
素 对 地 下 水 储量 的 影响 ,其 公式 为 : 
YCWSA, 
C = —— — —, i=1, +++, 12 (3) 


将 每 月 的 地 下 水 储量 变化 (GWSA ) 减 去 每 月 的 
气候 以 获得 地 下 水 储量 偏差 (GSD, Ground Water 
Storage Deviation) ,该 值 代表 地 下 水 储量 的 净 偏差 ， 
最 后 通过 减 去 地 下 水 储量 的 偏差 的 均值 ,在 除 以 其 
标准 差 得 到 GGDI。 


(GSD, E 


GGDI= (4) 
GSD 


式 中 :GSD, 是 在 GRACE 数据 提取 的 地 下 水 储量 变 
化 信号 中 剔除 月 气候 波动 得 到 的 地 下 水 储量 变化 
时 间 序 列 ; Xos 和 Seso 分 别 为 GSD, 的 平均 值 和 标 
准 差 。 
1.3.3 地 下 水 资源 量 可 持续 性 评价 及 指标 GGDI 
与 可 持续 性 指数 (Sustainability Index,SIT) 结 合 评估 
较 大 尺度 的 地 下 水 资源 量 可 持续 性 ,其 中 SI 为 地 下 
水 系统 可 靠 性 (Reliability, REL) 、 回 弹性 (Resil- 
ience,RES) 和 脆弱 性 (Vulnerability ,VUL) 的 函数 |. 
SI=RELxRESx(1- VUL) (5) 
根据 Thomas 的 分 级 标准 ,将 西北 干旱 区 地 下 水 
可 持续 性 分 为 5 个 等 级 : 低 可 持续 性 (0<SI< 0.2) , 较 
低 可 持续 性 (0.2<SI<0.3) ,中 可 持续 性 (0.3<SI< 
0.5) ,高 可 持续 性 (0.5<SI<0.75) , 极 高 可 持续 性 
(0.75<SI< 1), 
可 靠 性 (REL) 被 定义 为 系统 发 生 故 障 的 频率 。 
对 于 地 下 水 而 言 , 可靠 性 就 是 含水 层 储 存量 低 于 正 
常 条 件 的 历史 可 能 性 , 即 GGDI 处 于 可 靠 水 平 的 占 
比 。 将 西北 干旱 区 地 下 水 可 靠 性 分 为 5 个 等 级 : 低 
可 靠 性 (0<RELs0.25) , 较 低 可 靠 性 (0.25<RELs 
0.4) ,中 可 靠 性 (0.4<REL<0.6) ,高 可 靠 性 (0.6< 


REL«0.75) , 极 高 可 靠 性 (0.75<REL< 1). 

回 弹性 (RES ) 表 示 系 统 在 不 满意 的 情况 下 恢复 
到 满意 状态 的 速度 ,在 应 用 于 地 下 水 时 , 回 弹力 表 
示 从 不 满意 状态 恢复 到 满意 状态 的 可 能 性 , 即 
GGDI 从 负数 到 正 数 的 可 能 性 。 将 西北 干旱 区 地 
下 水 回 弹 分 为 5 个 等 级 : 低 回 弹性 (0<RES<0.2 ) , 
较 低 回 弹 性 (0.2<RES<0.3) ,中 回 弹 性 (0.3<RES<s 
0.5) , 高 回 弹性 (0.5<RES<s0.75) , 极 高 回 弹 性 
(0.75<RES< 1). 

脆弱 性 (VUL) 定 义 为 一 种 解释 不 满意 程度 和 出 
现 该 程度 的 概率 度量 。 对 于 区 域 地 下 水 系统 而 言 ， 
严重 程度 和 严重 程度 的 概率 即 地 下 水 干旱 指数 
GGDI 值 的 大 小 (Cs=GGDID) 及 其 出 现 的 概率 (e)。 可 
EN: 

VUL= se (6) 

将 西北 干旱 区 地 下 水 脆弱 性 分 为 5 个 等 级 : 低 
脆弱 性 (0<VUL< 0.1) , 较 低 脆弱 性 (0.1<VUL< 0.4)， 
中 脆弱 性 (0.4<VUL<0.6) ,高 脆弱 性 (0.6<VUL< 
0.75), 极 高 脆弱 性 (0.75<VUL< 1). 
1.3.4 Theil-Sen 趋势 分 析 与 Mann-Kendall 检 验 相 结 
合 ” Theil-Sen 趋 热 分 析 与 Mann-Kendall 显 著 性 检验 
是 2 种 非 参数 检验 方法 ,该 方法 不 要 求 时 间 序 列 必 
须 满 足 正 态 分 布 和 序列 自 相关 的 大 部 分 假设 ,而 且 
对 事件 序列 中 的 异常 值 不 敏感 ,能 够 有 效 处 理 小 的 
离 群 点 和 缺失 值 噪声 ,广泛 应 用 于 气象 和 水 文 的 时 
la Paap A, ASAI : 


X.,—X 
J 


p=wedil | Vj»i (7) 


I^ 
式 中 : B 是 所 有 数据 对 斜率 的 中 值 , B 的 正 负 表示 
时 间 序 列 中 的 趋势 方向 , 当 B>0 时 ,反应 上 升 趋 
势 ; 反 之 反应 下 降 趋 势 ; x 为 进行 分 析 的 指标 。 

本 研究 使 用 该 方法 , 取 置 信 水 平 w= 0.01 来 分 

析 地 下 水 储量 变化 (TWSA ) 的 时 空 变化 趋势 。 
1.3.5 最 大 信息 系数 MIC Reshef 等 在 2011 年 提 
出 了 最 大 信息 系数 (Maximal Information Coefficient, 
MIC) ,该 系数 可 用 于 挖掘 2 个 可 变 对 象 之 间 的 潜在 
关系 ,衡量 2 个 变量 之 间 的 线性 和 非 线性 关系 。 

其 原理 是 ,对 于 二 维 联合 随机 变量 (和 7) ,其 样 
ASE D-((x.y)i- 1,2,3, N) ,其 样本 容量 ; 
No WALD X A Y B5) (EL 33] 4 2 m I 个 不 同 的 
区 间 , 可 以 将 样本 空间 离散 化 为 mxn 的 网 格 G。 在 
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指定 网 格 下 可 以 进一步 估计 互信 息 1(Dl, )。 
I(DI,) - >, Ho pen y) | (8) 


p(x)pO) 
式 中 : DI, 表示 使 用 网 格 G 划 分 样本 集合 时 引入 概 
率 分 布 ;p(x),p0) 分 别 是 X 和 了 的 经 验 边 缘 密 度 p(x, 
是 和 了 的 经 验 联 合 概率 密度 。 
在 离散 化 样本 集合 D 时 ,同样 的 网 络 mxn 下 存 
在 多 种 不 同 值 域 划分 方式 ,所 有 可 能 的 网 格 C 上 的 
最 大 互信 息 记 为 


— 


I (D,m,n) 7 Max,I(Dl,) (9) 
进一步 标准 化 T (D,m,n) 
Min), e E eu) (10) 


in E log, min (m,n) 
基于 最 大 互信 息 , 随 机 变量 X,Y 之 间 的 MIC 定 
义 为 : 


Fr (a) 陆地 水 储量 变化 


TWSA/cm 


ED 


=f L 
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MIC(X) = Max, sy{M(X), )} (11) 
式 中 :BO 为 样本 个 数 的 函数 ,通常 设置 B(N)= 
N?* > MIC(X , Y ffc (ét vi ELO, 1] ,数值 越 大 线性 非 
线性 关系 相关 程度 越 大 。 


2 结果 与 分 析 


21 西北 干旱 区 地 下 水 储量 变化 

2.1.1 陆地 水 和 地 下 水 储量 变化 ”由 月 均 陆地 储 水 
量变 化 反 演 的 地 下 水 储量 变化 (图 2) 可 以 看 出 ， 
2002 一 04 一 2020 一 09, 西北 干旱 区 陆地 储 水 量变 化 
(TWSA ) 随 季节 变化 较 大 ,最 大 值 出 现在 2002 年 8 
月 ,其 等 效 水 柱 高 3.96 cm, 最 小 值 出 现在 2020 年 6 
月 ,其 等 效 水 柱 高 -4.42 em ,整体 呈现 出 周期 且 向 下 
的 趋势 。 地 下 水 储量 变化 (CWSA) 也 呈现 出 明显 


4 p (b) 地 下 水 储量 变化 
2L 
0r 


GWSA/cm 
b 
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图 2 陆地 水 储量 变化 和 地 下 水 储量 变化 


Fig. 2 Variations in of terrestrial water storage anomaly TWSA and groundwater storage anomaly GWSA 


的 .持续 的 下 降 趋 势 ,最 大 值 出 现在 2005 年 6 月 为 
3.06 cm ,在 2015 年 1 月 最 低 ,达到 -6.52 em 之 后 下 
降 趋势 有 所 减缓 ,开始 呈现 周期 波动 ,研究 期 间 等 
效 水 柱 高 累计 下 降 6.52 cmo 

21.2 地 下 水 储量 变化 验证 ”通过 对 研究 区 面积 
比 最 大 的 新 疆 地 区 进行 地 下 水 储量 变化 反 演 结 
的 验证 ,以 2003 一 2020 年 新 疆 水 资源 公报 的 地 下 水 
储量 年 鉴 数 据 作为 参照 (图 3) ,由 图 3 可 以 看 出 , 通 
it GRACE 反 演 获得 的 GCWSA 与 其 具有 相同 的 变化 
趋势 ,2014 年 新 疆 地 下 水 储量 最 低 ,为 443.9x108 
nm, 与 反 演 得 到 的 GCWSA 最 低 点 2015 年 1 月 相对 应 ， 
随后 在 2015 年 新 疆 多 地 发 生 连续 暴雨 洪水 灾害 ,地 


650- | —— GWSA -= 地 下 水 储量 -4 

eo. T^ | N E 
* [| i 
= a 10 g 
di 550 上 | | |] à 
E 500 上 = n X ] , © 
*) | Y s N J 

450 上 i 1-6 
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图 3 新 疆 地 下 水 储量 与 GWSA 
Fig. 3 Groundwater reserves in Xinjiang and 


groundwater storage anomaly GWSA 
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下 水 储量 有 所 上 升 *。 
2.1.3 地 下 水 储量 变化 空间 趋势 分 析 采用 Theil- 
Sen X% FI Mann- Kendall 显著 性 检验 相 结 合 的 方 
法 ,得 到 西北 干旱 区 地 下 水 储量 变化 趋势 (图 4)。 
研究 区 95.7% 区 域 在 2002 一 04 一 2020 一 09 呈现 出 
显著 变化 (as 0.01 ) ,其 中 约 有 19.3% 呈 现 出 显著 
的 正身 变化 , 即 地 下 水 储量 增加 ;80.7% 呈 现 出 显著 
负 向 变化 , 即 地 下 水 储量 减少 。 整 体 的 空间 趋势 呈 
现 出 3 层 条 带 式 分 布 , 自 北向 南 第 1 条 道 呈 现 增长 
趋势 ,增长 幅度 较 小 ,整个 条 带 东 西 跨度 最 小 ;第 2 
条 带 呈 现 出 明显 的 地 下 水 储量 减少 趋势 ,条 带 的 东 
西 跨度 最 大 ;第 3 条 带 呈 现 出 较为 明显 的 地 下 水 储 
量 增长 趋势 。 整 个 研究 区 大 部 分 区 域 地 下 水 储量 
变化 呈现 出 下 降 的 趋势 ,下 降 趋 势 最 显著 的 区 域 在 
Bur pod Hc .库尔勒 市 、 伊 宁县 交界 处 ,显著 增长 区 域 
最 明显 的 地 区 在 和 田 市 。 
2.2 地 下 水 干旱 指数 GGDI 

将 研究 区 基于 时 间 序 列 的 地 下 水 储量 变化 均 
值 ,去 除 基 于 GRACE 的 周期 变化 的 每 月 的 气候 值 


区 地 下 水 干旱 指数 GGDI>0, 最 高 值 出 现在 2003 年 9 
月 ,GGDI 为 1.95, 这 段 期 间 地 下 水 储量 较 长 时 间 处 
于 和 忽 余 状态 ,状态 较为 稳定 。2008 一 09 一 2012 一 10 
月 之 间 GGDI 在 0 上 下 波段 ,地 下 水 三 余 状态 和 赤字 
状态 相互 交替 。 在 2012 年 10 月 之 后 地 下 水 干旱 指 
数 GGDI 开 始 长 期 小 于 0, 且 呈现 减少 趋势 ,最 小 值 
在 2020 年 8 月 ,为 -1.66。 

2.3 地 下 水 资源 量 可 持续 性 评价 

23.1 可 靠 性 “可靠 性 是 西北 干旱 区 地 下 水 GGDI 
备 余 月 份 , 在 研究 期 间 总 月 份 中 的 占 比 情况 。 对 西 
北 干旱 区 地 下 水 可 靠 性 进行 量化 (图 6a) ,根据 可 靠 
性 分 级 标准 获得 西北 干旱 区 可 靠 性 分 级 结果 (图 
6b) ,以 及 西北 干旱 区 地 下 水 可 靠 性 各 等 级 占 比 情 
况 ( 表 1)。 

全 区 域 地 下 水 可 靠 性 平均 值 为 0.495, 达 到 中 可 
靠 性 水 平 。 其 中 ,中 可 靠 性 区 域 占 比 最 大 达到 
89.8%, 分 布 广泛 ; 较 低 可 靠 性 区 域 占 比 7.27% ,主要 
分 布 在 研究 区 东部 的 阿拉 善 左 旗 , 以 及 库尔勒 市 东 
南部 ,和 田 东 北部 ; 极 少 部 分 为 高 可 靠 性 区 域 ;无 低 


GRACE-C ( | 5b) ,得 到 地 下 水 储量 偏差 GSD( 图 
5c) ,再 将 地 下 水 储量 偏差 进行 归 一 化 后 得 到 基于 
GRACE 地 下 水 干旱 指数 GGDI( 图 5d)。GGDI>0 表 
示 地 下 水 储量 处 于 和 钥 余 状态 ,GGDI<0 表 示 地 下 水 
储量 处 于 赤字 状态 。 在 2008 年 9 月 之 前 ,西北 干旱 
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可 靠 性 区 和 极 高 可 靠 性 区 域 。 

2.3.2 回 弹 性 回 弹性 是 指 西北 干旱 区 地 下 水 
GGDI 从 赤字 状态 恢复 到 盘 余 状态 的 能 力 ,对 西北 
干旱 区 地 下 水 回 弹性 进行 量化 (图 7a) ,根据 分 级 标 
准 获得 西北 干旱 区 可 靠 性 分 级 结果 (图 7p) 以 及 回 
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4 利用 Theil-Senmedian 方 法 分 析 GWSA 时 间 变 化 
Fig. 4 Theil-Senmedian method was used to analyze the spatial variation of GWSA 


202206.00024v1 


chinaXiv 


ChinaXiv 合 作 期 刊 


10 T i T T d T f T T 7 T z 
(a) 地 下 水 储量 月 变化 | 


GWSA/cm 
© 


一 10 1 i 1 1 i 1 i 1 1 
2002-04 2004-11 2007—07 2010-02 2012-10 2015-05 2018-01 2020-08 


g 10 oT j ——T x 一 一 = 7—71 H 

S (b) 气候 周期 C | 

B o 

3... TT 
-10 


2002-04 2004-11 2007-07 2010-02 2012-10 2015-05 2018-01 2020-08 


10 一 一 一 一 一 一 — - -一 一 - - ~ -一 一 ——À] 


(c) 地 下 水 储量 月 偏差 | 
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图 5 基于 GEACE 的 地 下 水 干旱 指数 GGDI 
Fig. 5 Groundwater drought index GGDI based on GEACE 
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图 6 可 靠 性 评分 及 其 分 级 
Fig.6 Reliability score and its classification 
表 1 可 靠 性 等 级 占 比 弹性 各 等 级 占 比 情况 ( 表 2)。 
Tab.1 Proportion of reliability levels 全 区 域 地 下 水 回 弹性 平均 值 为 0.470 ,为 中 回 弹 
TREIA ul FTH 性 水 平 。 区 域 中 占 比 最 高 的 是 中 回 弹 性 区 域 ,达到 
低 可 靠 性 0<REL< 0.25 0 s ; " 
MAN 72.8% ,其 余 27.2% 为 高 回 弹性 区 域 , 主 要 分 布 在 西 
较 低 可 靠 性 0.25<REL< 0.40 7.27 z WAP RA 
ae pane in 北 干旱 区 南部 边界 沿线 ,以 及 西北 干旱 区 北部 ;中 
高 可 靠 性 0.6<REL<0.75 2.93 回 弹性 区 域 分 布 在 区 域 中 部 ,无 低 、 较 低 和 极 高 回 


极 高 可 靠 性 0.75<REL< 1 0 弹性 区 域 。 
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(b) 回 弹性 分 级 
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图 7 回 弹性 评分 及 其 分 级 


Fig.7 Resilience score and its classification 


R2 回 弹性 等 级 占 比 


Tab.2 Proportion of resilience levels 


回 弹性 分 级 范围 占 比 /% 
低 回 弹 性 0<RES< 0.20 0 
较 低 回 弹性 0.20<RES< 0.30 0 
中 回 弹 性 0.30<RES< 0.50 72.80 
高 回 弹 性 0.50<RES< 0.75 27.20 
极 高 回 弹 性 0.75<RES<1 0 


2.3.3 脆弱 性 ”脆弱 性 是 指 西北 干旱 区 地 下 水 赤字 
月 份 GGDI 的 程度 以 及 发 生 的 概率 对 西北 干旱 区 地 
下 水 可 靠 性 进行 量化 (图 8a) ,根据 分 级 标准 获得 西 
北 和 干旱 区 可 靠 性 分 级 结果 (图 8Bb) 以 及 脆弱 性 各 等 
级 占 比 情况 ( 表 3)。 

西北 干旱 区 平均 脆弱 性 为 0.404 ,脆弱 性 量化 结 
果 越 大 表示 区 域 越 脆弱 ,对 于 承受 极端 变化 的 能 力 
越 差 。 主 要 为 较 低 脆弱 性 和 中 脆弱 性 区 域 ,分 别 占 
53.69% 、42.2% , 其 余 小 部 分 为 高 脆弱 性 区 域 , 占 
4.1% ,脆弱 性 最 高 的 区 域 在 阿 瓦 提 县 和 库尔勒 交界 
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处 ,脆弱 性 水 平 达 到 0.67 ,无 低 脆弱 性 区 域 和 极 高 脆 
弱 性 区 域 。 
2.3.4 地 下 水 资源 量 可 持续 性 ”根据 获得 的 西北 干 
旱 区 地 下 水 可 靠 性 (REL), 回 弹性 (RES) 和 脆弱 性 
(VIL) ,获得 该 区 域 地 下 水 资源 量 可 持续 性 。 指 标 
中 的 最 小 值 对 该 指数 影响 最 大 , 当 3 个 指标 中 任 一 
值 为 0 时 ,区 域 不 具备 可 持续 性 。 对 西北 干旱 区 地 
下 水 资源 量 可 持续 性 进行 量化 (图 9a) ,根据 分 级 标 
准 获得 西北 干旱 区 可 靠 性 分 级 结果 (图 9b) 以 及 可 
持续 各 等 级 占 比 情况 ( 表 4)。 

西北 干旱 区 地 下 水 资源 量 可 持续 性 (SI) 的 全 区 
域 平 均值 为 0.28, 地 下 水 资源 量 可 持续 性 水 平 较 
低 。 区 域 总 体 地 下 水 资源 量 可 持续 性 总 体 呈 现 南 
北 两 侧 高 .中 部 低 的 态势 。 较 低 可 持续 性 区 占 比 最 
大 ,为 60.49% ,主要 分 布 区 域 东 南方 和 西北 方 , 包 括 
阿拉 善 左 旗 .哈密 市 .吐鲁番 市 等 ;其 次 是 中 可 持续 
性 区 占 比 为 32.47%, 主 要 分 布 在 区 域 西 南方 和 北方 
包括 和 田 市 、 阿 图 什 市 .喀什 市 .阿勒泰 市 等 ; 剩 下 
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图 8 脆弱 性 评分 及 其 分 级 


Fig.8 Vulnerability score and its classification 


202206.00024v1 


chinaXiv 


表 3 脆弱 性 等 级 占 比 


Tab. 3 Proportion of vulnerability levels 


脆弱 性 范围 占 比 /% 
低 脆弱 性 0<VUL< 0.10 0 
较 低 脆弱 性 0.10<VUL< 0.40 53.69 
中 脆弱 性 0.40< VUL< 0.60 42.20 
高 脆弱 性 0.60« VUL< 0.75 4.10 
极 高 脆弱 性 0.75< VUL< 1 0 


的 为 低 可 持续 性 区 占 比 , 为 7.03% ,大 部 分 分 布 在 阿 
瓦 提 县 ,零散 的 分 布 在 西北 干旱 区 中 部 ;无 高 可 持 
续 性 区 域 和 较 高 可 持续 性 区 域 。 
23.5 MIC 相关 性 评价 ”运用 python2.7, 基 于 最 大 
言 息 系 数 分 别 获 取 2002 一 2020 年间 西北 干旱 区 地 
下 水 资源 量 可 靠 性 、 脆 弱 性 、 回 弹性 和 可 持续 性 之 
间 的 最 大 信息 系数 MIC 以 及 皮尔 森 相 关系 数 。 
控制 地 下 水 资源 量 可 持续 性 的 直接 因素 是 地 
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下 水 储量 ,地 下 水 储量 变化 是 可 靠 性 、 回 弹性 、 脆 弱 
性 的 量化 依据 ,它们 分 别 描述 地 下 水 在 一 段 长 时 间 
序列 内 地 下 水 储量 低 于 正常 水 平 的 概率 、 从 低 于 正 
党 水平 恢复 到 正常 水 平 的 概率 、 低 于 正常 水 平 的 程 
度 与 其 发 生 概 率 的 乘积 。 

根据 相关 性 分 析 的 结果 表明 :可 靠 性 和 回 弹性 
与 可 持续 性 呈正 相关 ,脆弱 性 与 可 持续 性 呈 负 相关 
( 表 5)。3 个 指标 中 脆弱 性 与 可 持续 性 的 MIC 值 为 
0.512,R 值 为 -0.748,2 个 系数 的 绝对 值 均 大 于 另外 
2 个 指标 ,因此 脆弱 性 与 可 持续 性 相关 性 最 大 ,对 于 
地 下 水 资源 量 可 持续 性 的 影响 程度 最 大 。 


3 讨论 


综 上 所 述 ,2002 一 2020 年 西北 干旱 区 平均 地 下 
水 资源 量 可 持续 性 指数 SI 为 0.28, 可 持续 性 水 平 较 
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图 9 可 持续 性 性 评分 及 其 分 级 


Fig.9 Sustainability score and its classification 


表 4 可 持续 性 等 级 占 比 


Tab.4 Proportion of sustainability levels 


可 持续 性 范围 占 比 /% 
低 可 持续 性 0<SI< 0.20 7.03 
较 低 可 持续 性 0.20<SI< 0.30 60.49 
中 可 持续 性 0.30<SI< 0.50 32.47 
高 可 持续 性 0.50<SI< 0.75 0 
极 高 可 持续 性 0.75<SI< 1 0 


#5 相关 系数 分 析 


Tab. 5 Correlation coefficient analysis 


可 持续 性 (SIT) MIC R 

可 靠 性 (Rel) 0.235 0.084 
回 弹性 (Res ) 0.437 0.641 
脆弱 性 (Vul) 0.512 -0.748 


AK , BE Thomas 获取 的 塔里木 盆地 2002—2016 年 的 
地 下 水 资源 量 可 持续 性 指数 0.31 有 一 定 的 下 降 , 原 
因 是 2016 一 2020 年 以 来 ,新 疆 地 下 水 储量 持续 下 降 
(图 3),GGDI 长 期 处 于 赤字 状态 ,导致 地 下 水 资源 
量 的 回 弹性 和 可 靠 性 水 平 有 所 降低 ,进而 导致 地 
下 水 资源 可 持续 性 指数 的 降低 。 研 究 结果 与 吴 彬 
等 .中 基于 水 利 部 门 历次 水 资源 调查 评价 成 果 水利 
统计 资料 汇编 等 数据 分 析 得 到 的 新 疆 平原 地 区 
1956 一 2016 年 地 下 水 变化 趋势 一 致 。 

利用 GRACE 重力 卫星 结合 GLDAS 监测 地 下 水 
已 经 被 证 明 是 传统 地 下 水 监测 方法 的 有 效 奉 代 手 
段 ,尤其 是 在 基于 地 下 水 储量 变化 的 情况 下 ""。 不 
仅 降 低 了 传统 地 下 水 监测 作业 带 来 的 巨大 工作 量 ， 
还 克服 了 由 于 水 井 数量 少 分 布 不 均匀 造成 的 难以 
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进行 大 范围 长 时 间 序 列 监测 的 问题 。 相 较 于 传统 
的 GRACE 数据 ,研究 所 采用 的 JPL mascon RL06 数 
据 不 需要 进行 复杂 的 过 滤 和 其 他 后 处 理 , 约 束 性 较 
小 , 且 具 有 相似 的 精度 ,可 以 有 效 降 低 大 尺度 研究 
的 工作 量 。 研 究 时 间 跨 度 往往 是 影响 实验 结果 


人 类 活动 被 认为 是 地 下 水 资源 量 可 持续 性 变 
化 的 主要 驱动 因素 ,包括 地 下 水 开采 .地 下 水 灌溉 、 
城市 扩张 以 及 造林 工程 等 。 位 于 吐鲁番 西南 方 的 
托 克 逊 县 处 于 SI< 0.2 的 低 可 持续 性 区 ,受到 地 下 水 
开采 .农业 灌溉 影响 ,2002 一 2010 年 间 地 下 水 连续 


的 重要 因素 , 故 本 研究 选择 了 200 个 月 的 数据 ,时 间 
跨度 大 ,结果 可 信和 度 更 高 所 。 其 次 ,运用 基于 
GRACE 的 地 下 水 干旱 指数 GGDI 具 有 2 个 优点 :(1) 
GGCDI 能 直观 的 反应 地 下 水 的 赤字 和 盘 余 ; (2) 
GGDI 能 以 一 个 标准 化 数值 直接 反应 干旱 的 时 空 尺 
HE". Hate zer Hn C Bra rp RO. HB 
水 可 持续 指数 SI 综合 了 大 尺度 地 下 水 系统 的 可 靠 
性 , 回 弹 性 和 脆弱 性 来 评价 地 下 水 资源 量 可 持续 性 
其 结果 具有 较 高 的 参考 价值 ,可 以 为 地 区 地 下 水 保 
护 和 合理 利用 提供 参考 。 

影响 地 下 水 资源 量 可 持续 性 的 因素 主要 包括 
气候 变化 和 人 类 活动 。 降 雨 . 蒜 散 发 等 气候 因素 对 
地 下 水 储量 有 一 定 的 影响 。 比 较 2002 一 2020 年 , 西 
北 干 旱 区 年 内 与 年 际 降雨 与 GWSA 的 关系 (图 10)， 
在 年 内 尺度 上 由 GRACE 反 滨 的 地 下 水 储量 变化 与 
月 均 降 雨量 存在 明显 的 相关 性 ,在 降雨 量 最 大 的 7 
月 对 应 的 地 下 水 储量 变化 的 等 效 水 柱 高 最 高 。 在 
年 际 尺 度 上 , 2002—2020 年 ,降水 量 与 地 下 水 储量 
变化 的 等 效 水 柱 高 具有 相似 的 波动 ,如 2002 一 2008 
年 , 随 着 地 下 水 储量 变化 的 等 效 水 柱 高 随 着 年 降雨 
量 的 增加 而 上 升 , 随 之 降低 而 下 降 。 但 是 在 2002 一 
2020 年 ,地 下 水 储量 变化 的 等 效 水 柱 有 几 次 跳跃 式 
的 下 降 ,在 较 长 时 间 序 列 下 ,地 下 水 储量 变化 未 随 
着 年 降雨 量 的 变化 而 变化 ,说 明 还 存在 降水 量 以 外 
影响 地 下 水 储量 变化 的 驱动 因素 。 
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下 降幅 度 最 大 为 7.49 me 。 绿 洲 只 占据 了 西北 干旱 
区 不 到 10% 面 积 , 却 集中 了 干旱 区 95% 的 人 口 和 
90% 以 上 的 财富 2 ,绿洲 的 地 下 水 资源 相对 丰富 , 农 
业 几 乎 完全 依赖 地 下 水 灌溉 ,也 导致 绿洲 地 下 水 资 
源 量 可 持续 性 处 于 较 低 水 平 。 查 哈 尔 滩 绿洲 地 下 
水 资源 量 可 持续 性 处 于 0.2~0.3 的 较 低 可 持续 区 ,大 
规模 的 地 下 水 灌溉 ,导致 地 下 水 水 位 在 过 去 的 20 a 
平均 每 年 下 降 0.23 m^", Hao 4) HI TE BAR 
洲 和 敦煌 绿洲 的 研究 也 都 得 到 了 相似 的 结论 ,加 之 
近 十 几 年 以 来 西北 干旱 区 内 不 断 增 加 的 钻井 也 表 
明了 地 下 水 开采 、 农 业 灌 溉 是 导致 地 下 水 资源 量 可 
持续 性 降低 的 主要 因素 ;Wang 等 ”通过 主 成 分 分 析 
(PCA ) 比 较 河西 走廊 地 区 近 十 几 年 来 社会 经 济 组 合 

旧 标 与 气候 降水 指标 对 地 下 水 可 持续 性 的 影响 ,发 
现 前 者 解释 率 高 达 87.7% ,后 者 为 13.3% ,在 第 一 主 
成 分 中 人 口 因素 .CDP 等 与 地 下 水 可 持续 性 呈现 出 
负 相 关 性 ,但 影响 程度 较 小 ;西北 干旱 地 区 的 城市 
化 水 平 的 提升 也 是 地 下 水 可 持续 性 降低 的 因素 之 
一 ,根据 2000 年 .2010 年 .2020 年 Globeland 30 土 地 
利用 数据 统计 ( 表 6), 近 20 a 来 西北 干旱 区 人 造 地 
表 ( 即 建设 用 地 ) 面 积 增长 量 达 到 了 150% ,林地 面 
积 增长 量 为 40% ,耕地 面积 增长 量 21.5% ,城市 的 大 
规模 扩张 势必 造成 地 下 水 资源 量 可 持续 性 的 下 降 ， 
鲁 杨 “也 通过 夜间 灯光 数据 获取 的 西北 地 区 城市 
化 水 平 ,发 现 城市 化 水 平 与 地 下 水 储量 变化 存在 显 
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图 10 地 下 水 储量 变化 与 降雨 量 年 际 变化 .年 内 变化 


Fig. 10 Yearly and monthly variation of groundwater reserve and rainfall variation 
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表 6 2000 一 2020 年 西北 干旱 区 土地 利用 变化 
Tab.6 Land use change in arid area of Northwest China from 2000 to 2020 
土地 利用 类 型 /10'km 


耕地 林地 草地 灌木 地 湿地 水 体 人 造 地 表 裸 地 冰川 
2000 9.93 1.80 51.92 0.82 0.59 0.92 0.48 139.20 4.30 
2010 10.59 2.12 4141 1.50 0.55 1.11 0.56 148.02 4.11 
2020 12.06 2.49 41.03 1.33 0.98 1.37 1.20 145.84 3.67 


著 的 负 相 关 。 大 规模 的 造林 工程 也 是 影响 地 下 水 
资源 量 可 持续 的 重要 因素 ,2000 一 2018 年 毛乌素 沙 
地 植被 NDVI 显 著 提 升 ,地 下 水 储量 持续 下 降 与 ND- 
VI 呈 显 著 负 相关 ' 3, 另外 耕地 面积 的 持续 增加 也 势 


步 量 化 获取 了 西北 干旱 地 区 地 下 水 资源 量 的 可 靠 
性 、 回 弹 效 .脆弱 性 评价 分 析 西 北 干旱 区 地 下 水 资 
源 量 可 持续 性 指数 ,得 到 如 下 结论 : 

(1) 2002—04—2020—09 期 间 ,西北 干旱 区 地 


必 导 致 农业 地 下 水 灌溉 的 增加 ,也 是 造成 区 域 地 下 
水 资源 可 持续 性 的 下 降 的 因素 之 一 。 

对 于 西北 干旱 区 地 下 水 中 高 可 持续 水 平 区 域 ， 
以 节 流 保护 措施 为 主 :一 是 实行 地 下 水 取水 总 量 控 
制 与 水 位 控制 制度 ;二 是 明确 用 水 过 程 的 节约 要 
求 , 即 取 用 地 下 水 的 单位 和 个 人 遵守 取水 总 量 控制 
和 定额 管理 要 求 ,使 用 先进 的 节 水 技术 .工艺 和 设 
备 , 提 高 用 水 效率 ;三 是 强化 经 济 手段 的 运用 ,明确 
地 下 水 水 资源 税 费 的 征收 原则 ,四 是 细 化 地 下 水 保 
护 措施 ,防止 地 下 工程 建设 对 地 下 水 造成 重大 不 利 
影响 等 ;对 于 区 域内 较 低 、 低 可 持续 水 平 区 域 除了 
更 严格 的 开展 上 述 措 施 还 应 当 增加 水 资源 开源 措 
施 ,包括 开展 调 水 工程 等 ,遏制 地 下 水 资源 可 持续 
性 下 降 态 势 。 面 对 当前 西北 干旱 区 低 水 平 的 地 下 
水 资源 量 可 持续 性 形式 ,需要 加 快 落实 相关 规定 ， 
形成 一 个 合理 有 效 的 地 下 水 保护 管理 模式 ,促进 西 
北 干 旱 区 地 下 水 资源 可 持续 性 的 稳定 提升 。 

本 研究 局 限 性 在 于 没有 进一步 量化 考虑 影响 
西北 干旱 区 地 下 水 资源 量 可 持续 性 的 影响 因素 , 包 
括 上 自然 的 气候 降雨. 蒸 散 发 等 ,人 为 的 采矿 、 放牧、 
城市 扩张 、 防 护林 的 建设 等 ,是 下 一 阶段 主要 的 研 
究 方 向 。 男 外 ,由 于 地 下 水 资源 量 可 持续 性 评价 的 
结果 的 可 靠 性 . 回 弹性 脆弱 性 量化 时 的 主观 性 ,以 
及 数据 缺少 造成 不 同 版 本 数据 的 交替 使 用 ,使 研究 
结果 产生 可 接受 的 侦 差 ,存在 一 定局 限 性 。 


4 结论 


本 文通 过 GRACE 数据 和 GLDAS 水 文 模型 反 演 
西北 干旱 区 2002 一 04 一 2020 一 09 地 下 水 储量 变化 ， 
并 结合 基于 GRACE 的 地 下 水 干旱 指数 GGDI 进 一 


下 水 储量 整体 上 呈现 下 降 趋势 ,在 2015 年 1 月 等 效 
水 柱 高 最 低 ,达到 -6.52 cm 之 后 下 降 趋势 有 所 减缓 ， 
开始 呈现 周期 波动 ,研究 期 间 累 计 下 降水 柱 高 6.52 
cm。 空 间 上 看 ,区 域内 19.3% 的 区 域 地 下 水 储量 呈 
现 正 向 变化 ,80.7% 的 区 域 呈 现 负 向 变化 ,整体 自 北 
向 南 呈 现 条 带 式 分 布 ,中 部 条 带 地 下 水 储量 下 降 最 
为 明显 ,南部 条 带 地 下 水 储量 呈 上 升 趋势 。 

(2) 2002 一 04 一 2020 一 09, 西北 干旱 区 地 下 水 
干旱 指数 GGDI 变化 趋势 可 分 为 3 个 阶段 ,2002 一 
04 一 2008 一 08 GGDI 值 逐渐 减 小 ,但 值 大 于 0, 地 下 
水 储量 呈现 盘 余 状态 ;2008 一 08 一 2012 一 10 GGDI 
在 0 值 附 近 上 下 波动 ; 自 2012 年 10 月 开始 较 长 时 间 
处 于 赤字 状态 且 程 度 加 深 , 即 表明 西北 干旱 区 地 下 
水 干旱 程度 有 所 加 深 。 

(3) 通过 地 下 水 干旱 指数 GGDI 量 化 获得 西北 
干旱 区 地 下 水 的 可 靠 性 、 回 弹性 和 脆弱 性 以 此 评价 
西北 干旱 区 地 下 水 资源 量 可 持续 性 结果 显示 ,人 研究 
期 间 西 北 干旱 区 地 下 水 平均 可 靠 性 为 0.495, 回 弹性 
0.470, ,脆弱 性 为 0.404, 区 域 的 地 下 水 资源 量 可 持续 
性 指数 SI 为 0.28,3 个 指标 中 脆弱 性 与 可 持续 性 的 
最 大 信息 系数 最 大 为 0.512, 即 脆弱 性 对 于 地 下 水 资 
源 量 可 持续 性 影响 最 大 。 
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Evaluation of groundwater resource sustainability based on GRACE and 
GLDAS in arid region of Northwest China 


RUAN Yongjian", WU Xiugin™ 
(1. College of Soil and Water Conservation, Beijing Forestry University, Beijing 100083, China 2. Key Laboratory of 
State Forestry Administration on Soil and Water Conservation, Beijing Forestry University, Beijing 100083, China) 


Abstract: The arid region of Northwest China is a groundwater-dependent ecosystem. Groundwater resources are 
a strategic support for developing the arid region in Northwest China. It is important to understand the 
significance of groundwater changes in real time to reveal change characteristics of groundwater and regional 
economic development. Previously, traditional groundwater monitoring research was artificial and expensive and 
large-scale and long-term analyses were difficult. Research based on satellite information also lacks quantitative 
evaluations of groundwater sustainability. To explore recent groundwater sustainability in arid areas of Northwest 
China, GRACE gravity satellite and GLDAS data were used to invert groundwater storage changes and 
determine trend changes for 200 months from April 2002 to September 2020 in arid areas of Northwest China. 
The groundwater drought index, which is based on GRACE, was calculated to obtain the groundwater drought 
situation in the arid region of Northwest China. Finally, the GGDI was further combined with the sustainability 
index to quantify the reliability, resilience, and vulnerability of groundwater to evaluate the recent groundwater 
sustainability in the arid region of Northwest China. Results show that from 2002 to 2020, groundwater reserves 
were greatly reduced in the arid area of Northwest China, and the drought degree of groundwater has deepened. 
The average reliability, resilience, and vulnerability of groundwater is 0.495, 0.470, and 0.404, respectively, and 
the regional groundwater sustainability index (SI) is 0.28, which indicates a low groundwater sustainability level. 
This study reveals a temporal and spatial distribution and the evolution trend of groundwater sustainability in arid 
regions of Northwest China, which provides theoretical and supporting data for the protection and management 
of regional groundwater resources. 

Keywords: arid region of Northwest China; gravity satellite; groundwater drought index; groundwater resourc- 
es; sustainability 


